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We calculate and discuss the one-loop corrections to the photon sector of QED interacting to 
a background gravitational field. At high energies the fermion field can be taken as massless and 
the quantum terms can be obtained by integrating conformal anomaly. We present a covariant 
local expression for the corresponding effective action, similar to the one obtained earlier for the 
gravitational sector. At the moderate energies the quantum terms can be obtained through the 
heat-kernel method. In this way we derive the exact one-loop /3-function for the electric charge in 
the momentum subtraction scheme and explore both massless and large-mass limits. The relation 
between the two approaches is shown and the difference discussed in view of the possible applications 
to cosmology and astrophysics. 
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I. INTRODUCTION 

The interaction of electromagnetic field -with gravity 
is an important subject due to various astrophysical and 
cosmological applications. There are many publications 
on the subject and, in particular, many interesting -works 
devoted to the quantum corrections in the electromag- 
netic sector of the theory. In the early publications 
[1, 2], direct calculations -were performed through the 
Feynman diagrams and Sch-winger source method and 
later on, starting from [3], by using different modifica- 
tions of the heat-kernel approach [4-7]. The diagrams 
-were also used in the traditional style calculations [8, 9] 
in flat space-time and using a modified approach admit- 
ting a generalization to the curved space-time [11] [73]. A 
general consideration, including one and f-wo-loop calcu- 
lations and taking into account the temperature effects, 
has been given in [12]. 

One of the most important features of quantum cor- 
rections to the action of electromagnetic field is that they 
break the conformal invariance -which the theory pos- 
sesses at the classical level. Let us start by presenting 
a short list of the cosmological and astrophysical situ- 
ations where the violation of conformal symmetry can 
be relevant. (1) The breaking of conformal symmetry 
changes the equation of state of the radiation, that can 
be relevant in the radiation-dominated epoch in the early 
universe. The modified equation of state for the radiation 
enables one to construct interesting cosmological mod- 
els (see, e.g., [13] and references therein), including the 
ones [14] based on the quantum corrections of [3]. The 
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modified equation of state can produce the change in the 
la-w of expansion of the universe, entropy production and 
other phenomena. Also, this effect can slightly affect the 
red-shift dependence of both energy density of radiation 
and cosmic micro-wave background (CMB) temperature. 
Because of the gro-wing precision of astrophysical experi- 
ments, at some point this feature of quantum corrections 
can be relevant. (2) As a particular aspect of the pre- 
vious point, due to the broken conformal symmetry the 
rate of creation of the photons in the reheating period 
after inflation may be affected by quantum corrections, 
leading to the potentially observable consequences. (3) 
The violation of conformal symmetry is needed for the 
creation of initial seeds of magnetic field of the galax- 
ies, at the epoch of structure formation [15]. There is 
a -well-known attempt to explain this violation by quan- 
tum effects similar to conformal anomaly [16], during the 
inflationary epoch (see also [17]). It would be rather in- 
teresting to have a more complete understanding of the 
fleld theoretical mechanisms behind such violation, in- 
cluding the ones which can occur at much later epochs. 
The problem was further discussed, e.g., in [18, 19], and 
one can also see [20] and [21] for the recent reviews of 
the possible origin of the initial seeds of cosmic magnetic 
fields and related subjects. (4) Furthermore, the quan- 
tum corrections to the photon propagator can modify 
the position of the pole and hence produce the situation 
when the electromagnetic wave propagates with the ve- 
locity which is slightly different from the one in the purely 
classical case. A similar effect can take place in curved 
space and is sometimes characterized as a superluminal 
motion [22, 23] (see [24] for the recent review). In par- 
ticular, this effect may have a significant impact on the 
behavior of light in the vicinity of the black hole [25-27]. 
Indeed, a similar effect may also take place due to the 
presence of the boundaries and, in general, because of the 
macroscopic conditions [28] . Let us remark that the light 
propagation can be affected also by the nonlinear terms 
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[29] which can resuh from the quantum corrections and 
by the classical and/or quantum interactions to external 
gravity [30] or to other fields such as fc-essence [31]. It 
is obvious that the conformal symmetry forbids most of 
the possible changes in the wave equations and therefore 
the detailed study of the conformal symmetry breaking 
is relevant for this issue too. 

If thinking about conformal symmetry, the first re- 
markable observation is that in the realistic theories such 
as QED or the standard model (SM), and also in its gen- 
eralizations such as supersymmetric standard model or 
grand unified theories, the massless and conformal in- 
variant electromagnetic field couples to the other fields, 
which are all massive and therefore conformally nonin- 
variant. In the present paper we shall concentrate on the 
simplest case of QED and try to present the most gen- 
eral view of the violation of conformal symmetry within 
this theory. However, there is no qualitative difference 
other mentioned theories, which can be also considered 
elsewhere. 

One can formulate the two main questions: 

1) What is the mechanism of violation of the local con- 
formal symmetry in the intermediate energies and how 
one can link its violations at the ultraviolet (UV) and in- 
frared (IR) limit. Is it possible to derive such relations? 
This is not a purely technical question, because during 
its evolution the Universe is passing different phases and 
it is desirable to know how the conformal symmetry is 
violated not only at the inflation epoch or at the present 
epoch, but also around the recombination epoch or at 
the time period when the cosmic structure starts to form. 
The main purpose of the present paper is to explore this 
issue. 

2) To which extent the finite quantum corrections and, 
in particular, the violation of local conformal symmetry, 
are universal? In other words, do we have some ambigu- 
ity in the quantum terms? Answering the last question is 
of course very significant, because the effective action, in 
general, is not a uniquely defined object. Usually, it de- 
pends on the choice of parametrization of quantum fields, 
in particular on the gauge fixing choice, on the calcula- 
tional schemes, regular ization, renormalization, etc. 

Our purpose is to derive the most general expression 
for the one-loop quantum correction to the electromag- 
netic sector of QED in curved space-time. As we have 
already mentioned above, the physical situations for such 
quantum corrections in the early universe and at the later 
period are very different. In the first case the fermions 
can be treated as almost massless. In this case the mech- 
anism for the violation of conformal symmetry is the well- 
known conformal anomaly. The advantage of conformal 
anomaly as a method of deriving quantum corrections 
to the classical effects is its simplicity, direct relation to 
the UV divergencies and consequent universality. At the 
opposite end of the energy scale the masses of the vir- 
tual fermions are much greater than the energies of the 
real photons or, equivalently, the energies of the external 
tails in the loops of massive fermions. In this case we ob- 



serve the decoupling of the quantum contributions of the 
massive fields, according to the Appelquist and Caraz- 
zone theorem [32] . The violation of conformal symmetry 
still exists, but it is related to the remnant higher deriva- 
tive terms in the effective action, which are quadratically 
suppressed by the fermion masses. Also, in curved space- 
time, there may be curvature-dependent terms which can 
break down the local conformal symmetry. Since we in- 
tend to study the curved-space theory, we will always 
concentrate our attention on the local conformal symme- 
try. 

The paper is organized as follows. In the next sec- 
tion we consider the simplified version of QED with a 
massless fermion and obtain the noncovariant and also 
nonlocal covariant forms of the anomaly-induced effective 
action. Some part of this section was previously known, 
but we present these results for the sake of completeness, 
and also to establish the most general local representa- 
tion for the covariant expression. In Sect. 3 we start to 
consider massive case and calculate the first three coeffi- 
cients of the Schwinger-DeWitt expansion in an arbitrary 
space-time dimension D. Even though the focus of our 
interest is mainly on the AD case, it proves useful to have 
a more general _D-dimensional result. The calculations 
are performed in two different schemes, such that we can 
try the limits of universality of the quantum corrections. 
The scheme dependence of the one-loop effective action 
is explored in a parallel paper [33]. In Sect. 4 we use the 
heat-kernel method to calculate the complete one-loop 
formfactors for the electromagnetic field sector. In the 
same section we consider the UV limit and establish the 
relation of the massless limit for the one-loop corrections 
for the massive case with the ones derived via confor- 
mal anomaly. In Sect. 5 we consider the renormalization 
group equations for the massive case and derive different 
forms of the low-energy decoupling law corresponding to 
the distinct calculational schemes. In Sect. 6 we an- 
alyze the corresponding running of the effective charge 
and compare it to the one in the minimal-subtraction 
renormalization scheme. Finally, in the last section we 
present some discussions and draw our conclusions. 

II. ANOMALY-INDUCED ACTION FOR THE 
METRIC AND ELECTROMAGNETIC FIELDS 
BACKGROUND 

Let us start from a brief survey of massless confor- 
mal QED. The theory should be formulated in curved 
space-time and hence the action depends on the electro- 
magnetic potential A^, Dirac spinor field V' and on the 
external metric g^^- As far as we are interested in the 
local conformal symmetry, one of the useful parametriza- 
tions of the metric is 

9tiv = 9iiv -e^" , a = (j{x) , (1) 

where fjf^u is the fiducial metric with fixed determinant. 
For example, in the case of the cosmological metric, using 
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spherical coordinates, we have 
g^, = diag (l, - YZ^' 
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Separating a{x) in (1) proves to be a useful tool, espe- 
cially because of the relation 
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which is valid for any functional of the metric 

and maybe other fields. If we replace the action of some 
theory in curved space at the place of ^[3^1^], then the 
l.h.s. of the above relation is nothing else but the trace 
of the corresponding energy- momentum tensor, Tj^. In 
order to remove the effect of other field variables, it is 
sufficient to use the corresponding equations of motion. 

The vanishing trace of the energy-momentum tensor 
implies that the conformal factor of the metric decouples 
from the matter. However, the situation changes dramat- 
ically if we take quantum effects onto account. In the last 
case the corresponding theoretical phenomenon is called 
the trace anomaly [34] (see also [35, 36] for review and 
many further references). 

The classical action of electromagnetic field is 



(3) 



and possesses local conformal invariance. The last means 
this action docs not change under simultaneous transfor- 
mation of the metric and of the vector A^, namely 
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Let us note that the difference between conformal weight 
and dimension for the vector field is due to the vector 
field definition in curved space-time, 

A^ = A,el, e>^7?afc=.9^., el g^-" = rj'^" . (5) 

We are interested in the corrections to the action (3) 
due to quantum effects of the fermion 

Sf = i j {V'7''(V^-*eA^)^ - imi^V} -(6) 

The conformal transformation rule for spinors is 

The metric is always transformed like in (4). Indeed the 
action (6) is conformal invariant only when the spinor 
mass is zero, m = 0. All those fermions which couple to 
the electromagnetic field as in (6) are massive, however, 
the relevance of the mass terms depends on the energy 
scale. For instance, if we are interested in the quantum 
effects of fermions close to the infiationary epoch, the ki- 
netic energy of the real fermions and (more important) of 
photons is much greater than the mass of the fermions. 



In this situation taking a massless spinor field is a legit- 
imate approximation, so let us start from this case and 
take m = 0. 

The derivation of conformal anomaly in the presence 
of background metric and electromagnetic field has been 
discussed before [34, 35], and we can use the known re- 
sult. The conformal anomaly can be used to construct 
the equation for the finite part of the one-loop correc- 
tion to the effective action of the background metric and 
electromagnetic potential, 
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where 



— — ^ixva^ ~ 2i?^^ -|- (1/3) i?^ 

is the square of the Weyl tensor and 
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is the integrand of the Gauss-Bonnet topological term. 
The coefficients uj, b, c depend on the number of scalar 
Ns, fermion Nf and massless vector A^^ fields as 

111 

w = — Ns^ — Nf^ — N^, 

120 ^ 20 ^10 
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Ns Nf iV^ , 

360 360 ' 180 

180 30 ^ 10 



(8) 



Also, ^ depends on the number of charged scalars (in 
case of scalar QED) and spinors. As far as the massless 
approximation can be applied in the very early universe, 
the number of fields which contribute to these coefficients 
is not necessarily restricted by the QED framework. 

The solution of Eq. (7) is straightforward [37] (see 
also generalizations for the theory with torsion [38] and 
with a scalar field [39]). The simplest possibility is to 
parametrize metric as in (1), separating the conformal 
factor a{x) and rewrite Eq. (7) using (2). The solution 
for the effective action is 



r — Sc\gfiv, Aij\ 



1 



(47r)2 

+ 0aFl^ + ba{E - |ni?) + 260-A40- 
- ^(c+^6)[i?-6(W)2-(na)]2} (9) 

where S'c[.9^;y, A^,] = Sc[gfiv] is an unknown conformal 
invariant functional of the metric and A^^ , which serves 
as an integration constant for the Eq. (7). All quanti- 
ties with bars are constructed using the metric g^i,, in 
particular 

^ ixv ~ ^ nvJ^ajsg g 
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Furthermore, A4 is the fourth derivative conformally co- 
variant operator acting on dimensionless scalar 

A4 = + 2R^'"V^Vu - ^ iin + ^ i?;^ . (10) 

The solution (9) has the merit of being simple, but an 
important disadvantage is that it is not covariant or, in 
other words, it is not expressed in terms of an original 
metric g^v In order to obtain the nonlocal covariant 
solution and after represent it in the local form using 
auxiliary fields, we shall follow [37, 40]. The presence 
of the F'^j, terms does not require any essential changes 
compared to the consideration presented in [36], in par- 
ticular this term can be always taken together with the 
one. So, we present just the final result in the nonlocal 
form, which is expressed in terms of the Green function 



G{x,y) of the operator (10), 

^4,xG{x,y) = 6{x,y) . 

Using the last formulas and (2) we find, for any 
Mdniy) = ^iQiJ-u e^'^) , the relation 

= 4VIA4 A = 4V5A4 A . 
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In particular, we obtain 

^induced ~ ~t~ Tfo ~t" j 

where 



(11) 



(12) 



= J 1 d^x./^ I dS^/^){ujC^ + pFl,)^G{x,y) (e-^UB)^, (13) 
^b = \ j d^x^^) j dS^){E-'^-UR)^G{x,y){E-'^UR)^ (14) 

and 



The nonlocal expressions for the anomaly- induced ef- 
fective action can be presented in a local form using two 
auxiliary scalar fields ip and ip [40]. Let us give just a fi- 
nal result which has an extra electromagnetic terms com- 
pared to the one described in [36] 



r = 5e[5M., ^a] 
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(16) 



The local covariant form (16) is dynamically equivalent to 
the nonlocal covariant form (12). The complete definition 
of the Cauchy problem in the theory with the nonlocal 
action requires defining the boundary conditions for the 
Green functions G{x,y), which show up independently 
in the two terms (13) and (14). The same can be 
achieved, in the local version, by imposing the boundary 
conditions on the two auxiliary fields <f and ijj- 

Let us separate the part of the effective action (16), 
which has direct relation to the electromagnetic terms 



and therefore represents a one-loop correction for the 
classical action (3), 



r = Sc[g^^,Ax] + 



d^x 



Stt 
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(17) 



Let us note that the presence of the Ay^-independcnt term 
'i^[E — (2/3) Di?] is relevant, because only this term pro- 
vides violation of local conformal symmetry in the whole 
expression. For instance, the connection with the nonco- 
variant presentation (9) is through the relations 



and 



s/g(E-'^OR 



^{e-'^OR + AKag). (18) 



Another important observation is that the anomalous 
metric dependence of the F^^^ in Eq. (17) appears due 
to the coupling of F'^^ with the auxiliary scalars ip and 
ip. This dependence is nontrivial because these two fields 
have different space-time behavior due to their distinct 
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dynamical equations and independent initial and bound- 
ary conditions. 

As an example of the use of the auxiliary fields ip and 

(p, we mention that the different choices of their boundary 
conditions enable one to classify the vacuum states for the 
semiclassical black hole [41]. It would be interesting to 
explore the initial and boundary conditions for ip and tp 
such that it could be applied, for example, for calculating 
the creation of the seeds of a magnetic field during infla- 
tion. However this issue requires (and deserves) a special 
detailed investigation which goes beyond the scope of the 
present paper. 



III. MASSIVE FERMION CASE: 
SCHWINGER-DEWITT EXPANSION 

The results of the previous section can be seen as fol- 
lows: in the case of massless and conformal fields we can 
calculate an important (sometimes the most important) 
part of the one- loop effective action by using conformal 
anomaly. The qualitative explanation of this possibility 
is that the conformal anomaly is essentially controlled by 
the logarithmic divergences of the theory. And, on the 
other hand, the logarithmic divergences are intimately 
related to the UV behavior of the theory, or, better to 
say, to its reaction to the scaling in the high energy re- 
gion. In the massless theory there is no natural scale and, 
therefore, all kind of quantum processes can be actually 
seen as "high energy". For this reason the UV diver- 
gences provide much more information on the finite part 
of effective action in the massless conformal case than 
they do for massive fields. The question is what we can 
do in this case. 

The theory of a massive quantum field has the nat- 
ural scale established by the mass. In this situation to 
notions of "high energy region" (UV) and "low energy re- 
gion" (IR) assume some relations between the energy of 
the field and the mass of the quantum field. The quan- 
tum effects of massive fields are supposed to be close 
to the ones of the massless fields in the UV and follow 
the decoupling theorem [32] in the IR. In both cases the 
most relevant finite part of effective action is given by 
nonlocal expressions, but the structure of nonlocalities 
are rather different in the two cases. As far as we are 
interested in the photon propagation, all that we need 
is a formfactor of the electromagnetic term. We shall 
present this result, at the one-loop level, in the next sec- 
tion. In the present section we consider the local terms 
which correspond to the Schwinger-DeWitt expansion of 
the effective action in the massive case. The analysis of 
the coefficients of this expansion will prove useful for bet- 
ter understanding of the structure of the formfactors of 
our interest, moreover deriving these coefficients provides 
an efficient check of correctness of the consequent general 
calculations. 

The one-loop effective action (EA) in the metric and 
electromagnetic sectors can be defined via the path inte- 



gral 

^r[g^.,A^] ^ J jj^jj^^iSQED ^ (19) 

where 

Sqed = Sf + Sem (20) 

and the actions Sf, Sem are defined in (3) and (6). Since 
the action Sqed is bilinear in the spinor fields, we find 
(see, e.g., [42]) 

= - 1 Ln Det ff, (21) 

where 

H = ^(^^^v^-^m-ie^^'A^) (22) 

is the bilinear form of the action (6). Here and below we 

assume but usually do not write explicitly the identity 
matrix 1 in the space of Dirac spinors. 

The derivation of (21) can be performed by several 
methods. Here we intend to use the heat-kernel ap- 
proach and the Schwinger-DeWitt technique. For this 
end we need to reduce the problem to the derivation of 
Ln Det O, where the operator O should have the form 

6=0 + 2/i''V^ + H . (23) 

An obvious way to achieve the desired form of the opera- 
tor is to multiply H by an appropriate conjugate operator 
H*, 

6 = H -H* (24) 

and use the relation 

Ln Det H = Ln Det 6 - Ln Det H* . (25) 

It is clear that one can obtain the desirable form of the 
product by simply choosing H* = H. In this case we can 
obtain the result (21) by taking 

Ln Det H = ^ Ln Det . (26) 

At the same time there arc many other possible choices. 
For instance, the calculations can be performed in the 
most simple and economic way if we choose the conju- 
gated operator with the opposite sign of the mass term, 

= _ i (^/^v^ + - 167"^^) . (27) 

According to [43] , the contributions of the two operators 
are identical, 

Ln Det H = Ln Det .^i , (28) 
such that we can still use the formula similar to (26), 

Ln Bet H = ^ Ln Det (HHi) . (29) 
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The expressions for the elements of the operator (23) in 
this case are rather simple, namely 

III = -\R + m^-'-^ri''F^. 

- ie{V^A^) - ie^A^'A^. (30) 

Another possible choice of H* is 

4* = -^(7''V^-^m). (31) 
In this case we have to use the relation 

Ln Det ^ = Ln Det {HHl) - Ln Det . (32) 

It is obvious that the contribution of Hi does not depend 
on Afi, just because this operator docs not depend on ^4^. 
Therefore, for evaluating the contribution of Ln Det H 
we can take the A^-dependcnt terms from the contri- 
bution of the product [HHf) with the coefficient one, 
while for the ^^-independent terms the coefficient must 
be one-half. 

The UV divergences derived within the two schemes, 
the first one with H* and the second one with are 

equal to each other. At the same time there is an essen- 
tial difference between the finite parts of the correspond- 
ing effective actions derived within the two schemes. The 
general discussion of this difference (called multiplicative 
anomaly, [44, 45]), is discussed in the parallel paper [33]. 
Let us present the technical part of the consideration, 
which are also relevant for our main targets, that is vio- 
lation of conformal symmetry and, in general, quantum 
corrections. 

The elements of the operator O for the case of the 
operator H2 have the form 

% = ~R + m^ + eM^''A^. (33) 

It proves useful to calculate the first Schwinger-DeWitt 
coefficients in an arbitrary space-time dimension d. For 
this end we will need the expressions for the operators 

P = n + ^ i? - V^/i" - ft^/i" (34) 

and 

S1J.U = R^n - V^/i^ + V^h^ - h^h„ + Khfj, , (35) 

where R^v = l\W^,Wu\ is the commutator of the two 

covariant derivatives acting on the corresponding fields. 
In our case of Dirac spinors space (sec, e.g., [42]), 

R-Hu V' = [V^ , V,] V = J RtJ^i^Xr 7^7^ , 



7^7" V^V. = ^-\r. (36) 

It is important to remember that formulas (34) and (35) 
do not depend on the dimension d, and the general ex- 
pressions for the coincidence limits of the Schwinger- 
DeWitt coefficients do not depend on d either. 

The expressions for P and S^j^v for the two calcula- 
tional schemes (with and correspondingly) are 
as follows: 

Pi = m'-^R-'-^^'^YF,^, (37) 
Si,ij,u = RoiHixv + ieFij,^; 

versus 

A = - ^ ^ - 7 ^^t" F^,, + emYA^ 

+ '-liWA,)-'^^^A'^A^, (38) 

S2,i^t. = -^7"7'^-Ra/3/i.. + ^7"(7^V^A„-7^V^Aa) 

g2 

+ ^ 7a (7^7/37.^ - 7-/7/37/^) ^"^'^ • 

The operators (37) and (38) enable one to calcu- 
late the coincidence limits of the first coefficients of the 
Schwinger-DeWitt expansion, by just using the known 
general expressions (see, e.g., [46, 47]). 

Let us start from the first nontrivial coefficient. For 
the sake of simplicity we shall use the notation 

ttfe = Tr lim ak{x,x') 

x'—^x 

for the functional traces (including space-time integra- 
tions) of the coincidence limits. We know [46, 48] that 
for the operators of the form (23), the first coefficient is 
given by the expression 

ai = sTrP = ~ j d'^Xy/^ tr P , 

where the operators P arc given by the expressions (37) 
and (38). The symbol sTr means the functional trace, 
which is taken according to Grassmann parity of the cor- 
responding functional matrices. Some pedagogical ex- 
amples, including the operators with both bosonic and 
fermionic blocks can be found in [42]. For the general 
rf-dimensional case we obtain the expressions 

ai^i = - J d^'x^iAm^ -^Ry (39) 

and 

ai,2 = - J d'^x^{Am^-^R + 2ie{V^A^) 

- {d-2)e^A^'A^Y (40) 
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The first expression (39) is independent of the electro- 
magnetic potential and is, therefore, gauge invariant. 
However, the second expression is obviously different in 
both respects, in particular it becomes gauge invariant 
only in the d = 2 case (in what follows we shall call it 2d, 
and exactly the same way for other dimensions, e.g. Ad 
means d = 4 etc). 

How can we understand the difference between the re- 
sults coming from the two calculational schemes, based 
on and i?2 conjugate operators? Let us note that the 
operation of introducing the conjugate operator (24) can 
be seen in two different ways. 

First, we can see it as a change of the quantum vari- 
ables in the path integral (19), that is taking ip = H*x, 
where x is a new quantum variable. In this case we have 
to take into account the contribution of the functional Ja- 
cobian of such transformation. This Jacobian is in fact 
related to another path integral, similar to (19), namely, 
for the second case, 

J DxDx exp [i J d^x^xH^x} ■ (41) 

This integral does not depend on the electromagnetic po- 
tential, but only on the metric, and hence it is gauge in- 
variant. However, the result of the change of variables is 
quite different, for one meets another functional integral, 

J DxDij e^p{ij d^x^ i> {H ■ H;) x} • (42) 

This integral does not possess gauge invariancc and hence 
it is not a surprise that the result of the change of vari- 
ables also does not possess it. 

Another possible understanding of the operation (24) 
is through the well-known relation 

Ln Det {H* ■ H^) = Ln Det H* + Ln Det H2 ■ (43) 

The two terms in the r.h.s. arc gauge invariant due to 
the reasons explained above, namely the first term is a 
result of the gauge-invariant functional integration and 
the second term simply does not transform, just because 
it does not depend neither on neither on spinor field. 
Hence if we find the violation of gauge symmetry in the 
l.h.s. (as we actually did), this indicates the violation of 
the "identity" (43). 

The relation similar to (43), 

det {A -3)= det A - detB 

can be easily proved for the finite-size square matrices 
A and B. However, the proof can-not be generalized 
for the differential operators, which have an infinite-size 
matrix representations. In fact, some mathematicians 
and physicists were looking, for a long time, for an ex- 
ample when this relation does not hold [44, 45], by us- 
ing the C-regularization technique [49] . The phenomenon 
of such possible violation has been called Multiplicative 
Anomaly. However, the results of the mentioned works 



met justified criticism [50-52] because it is in fact difficult 
to distinguish the effect from the usual renormalization 
ambiguity. 

From the perspective of Multiplicative Anomaly we 
can observe, that the key relation (43) hold in 2d and 
only in 2d. However, this dimension is very special for 
the ai, because in this particular dimension the coeffi- 
cient ai defines the logarithmic divergence of the theory. 
At this point we conclude that, in the case under consid- 
eration, the relation (43) does hold for the logarithmically 
divergent part, but may be violated in the other sectors 
of the effective action. In a moment we shall see that this 
statement is also valid for the next two coefficients. 

The general expression for the second Schwinger- 
DeWitt term (called also the "magic" coefficient) is 

02 = sTrd2{x,x) = sTr{^{Rl^„f^-Rl^ + nR) 

+ IP' + li^'P^ + Y^^l-^}- (44) 

Direct calculations using formula (44) and the expres- 
sions (37) and (38) yield the following results 

02, 1 = ^ j d'^x^ {ASe^F^^F''-' + R^ - 2ARm^ 
- SR^^a^R^""!^ + lUm''). (45) 

for the first calculational scheme based on and 

a2,2 = ^ / d''x^[2Ade^{V^An{yM 

- 96 e'^(yA''){V^Af,) + R^ + 144 m"* - 24:Rm^ 

- 3 R^^affR'"'"^ + 12(d - 4)(i? - 12m^)A^'A^ 

+ 6{d-2){d-4)A^A^'A^A•'e^^ (46) 

for an alternative calculational scheme based on ifg. It 
is easy to see that the two expressions for a2 presented 
above follow the same pattern as the ones for ai which 
we considered earlier. Namely, the expression for a2, 1 is 
gauge invariant independent of the space-time dimension 
d and, in fact, its dependence on d is rather simple. On 
the contrary, 02,2 manifest much more complicated de- 
pendence on d. This expression is gauge invariant only 
in the M case, when it also coincides with 02,1. We 
know that this coincidence is because the 02 defines the 
logarithmic divergences in Ad and only in Ad. So, we 
observe that the logarithmic divergences of the theory 
are scheme independent and that for the divergent sector 
of the theory, in dimensional regularization, the relation 
(43) actually takes place. At the same time, for any other 
dimension d 7^ 4, the two expressions are dramatically 
different. We can indicate two aspects of this difference. 
First, the 02, 2 (^ 7^ 4) expression is not gauge invariant. 
This means the quantum contributions depend not only 
on the transverse part of A^, but also on the longitudinal 
part of this vector. Let us note that the 02, 1 (rf ^ 4) is 
perfectly gauge invariant and no longitudinal propaga- 
tion takes place in this case. 



8 



Finally, as a final test of the difference between the 
two schemes of deriving the photon propagator, let us 
calculate the 03^ 1 and 03, 2 coefficients. For this end we 
will need the general expression for the 03 coefficient, 
which has been derived by Gilkey [53], and checked by 
Avramidi [54] who also derived 04 coefiicient. This gen- 
eral expression is rather bulky and we simply refer an 
interested reader, e.g., to the Eq. (2.160) of [54]. One 
observation is in order. The result of [53, 54] correspond 
to the operator of the more simple form 

O = Q + n . (47) 

if compared to our (23). However the way to generalize 
from (47) to (23) is already known for a long time from 
[47], where it was applied for the 02 coefficient. One has 
to start by constructing the generalized covariant deriva- 
tive acting in the fermionic space 

It is easy to see that the operator (23) can be presented 
in the form (47) with the new derivative and, also, II 



replaced by P. The last thing to do is to replace the 
commutator R^i, by S^j^^, defined in (35). Then the 
Schwinger-DeWitt technique can be developed with the 

covariant derivative instead of V^, and eventually 

gives the result in terms of P, Sf^i, and curvature tensor. 

Consider the most simple ff^-based calculational 
scheme. Using the original approach for the operator 
(47) we arrive at the already known result of Ref. [3], 

03,1 = ^(2i?^,„/3F'^''i^"'3 - 26i?„,F'^^F^" 
+ 24 V^F^" VaF^" + 5R F^"" F^, ) . (48) 

It is easy to see that this expression is gauge invariant 
independent of the space-time dimension d. Let us now 
see what is the situation with the second choice of the 
calculational scheme, based on the operator In this 
case one has to use the approach described above, that is 
to deal with the operator (23) and use generalized covari- 
ant derivative D^. After an involved algebra we arrive 
at the expression 



03,2 



AA " 2880 r^°^^'^^°^^^^~^°^''"°^''"~^^^"^''"'^"''^''" + 24(n^")[(rf-3)(nA„) + 2V„(VA)] 

- 24(V„V^A^) [(V'V^A") - (V^V^A'^)] + [(18 - 7d)i?2^„^ - 8(9 - d)Rl, - 6(5 - d)R^] 



+ 8Rnua0 [4(V"A'')(V''A'^) - 8F^"'F"'^ - 3{d - 4)(V''^")(V''^^) - R^'^'^'^A^Ax + lOR'''^ A"" A"] 
+ 16R^^[10{WA){V^'A'') + (V"A'')(5VaA'' - 2V"A„) - {d - 5){V^'A"){V'A^) - 2R^'^A"A-'] 
+ WR[2{d - 5)iVMi^'"^l - 2(VA)2 + SF^.F'"' + 2R^,A''A^] - 24(^^7?) [(V.A«A„) - (V^A^A,)] 

(49) 

I 



12{d - 2)A°'Ac,aR - 48(V"i?^^.a/3)(V''A'^^'')} . 



The last expression is rather cumbersome and difficult 
to deal with. Apparently it is not gauge invariant and 
is different from (48). In order to complete the analysis, 
we need to answer two questions: (i) How we can prove 
that the two expressions are distinct or identical in one 
or another particular dimension? (ii) Is it true that 
the two expressions (48) and (49) are distinct, in one 
or another particular dimension? In case of as we are 
mainly interested in the 6d case, where one can hope to 
find the gauge invariance of (49) and also the equivalence 
between (48) and (49). 

Let us start by making the most simple test. Con- 
sider the propagation of the purely transverse part of 
Afj,, that means to collect the terms which contain the 
terms A^d^A^ in both expressions. It is important that 
the transverse terms can-not be affected by the possible 
violation of gauge invariance in the expression (49), so 



this particular test is independent from of rest. The rela- 
tively simple calculation shows that the two propagators 
are identical in 6d and only in 6d. This solves one-half 
of the item ( ii), namely we learn that out of 6d the two 
expressions (48) and (49) are different. 

We could not precisely prove that the two expressions 
are identical in 6d for an arbitrary metric background. It 
is likely, however, that this is the case. The arguments 
in favor of this conclusion come from the analysis of a 
particidar simple metrics. Let us consider the case of 
de Sitter space-time, where one can actually check the 
equality of Eqs. (48) and (49). Even in this case it is 
not easy to work with (49), but one can perform some 
qualitative analysis of the expression for aa, 2- First of all, 
by direct inspection we can verify that the terms without 
scalar curvature are the same in both expressions. 

As a next step, consider the terms proportional to the 
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scalar curvature, which can be symbolically identified as 
RFF. In order to analyze these terms, one has to take 
into account also the terms proportional to the deriva- 
tives of the curvature tensor components in (49). For 
definiteness, we call them DR-terms. If one integrates 
these terms by parts, it turns out that some of them do 
contribute to the i?^j4^-structures. Since these terms can 
be written as total derivative terms in de Sitter space- 

(2) 

time, we can introduce them in Og with an arbitrary 
coefHcient. To make the correct choice of this coefficient, 
one has to additionally evaluate those terms which are 
proportional to B? and choose the coefficient of the 
DR-terms such that they do cancel the E?A^ terms in 
(49). After that, the terms proportional to RFF, gener- 
ated by integration by parts of the DR-terms, sum with 
the ones which come from (49) and finally turn out to be 
the same as (48). We do not include the corresponding 
formulas here, because they are rather boring. 

Finally, we have a strong argument to believe that 
qualitatively the situation with 03 is the same as with 
the ai and 02 coefficients. All these coefficients are uni- 
versal (scheme independent) in the dimensions 2, 4 and 6, 
correspondingly. And, at the same time, they are essen- 
tially scheme dependent in other dimensions. It is easy to 
see that this fact implies that, in any dimension, the log- 
arithmic divergences are scheme independent and gauge- 
invariant. At the same time, quartic and quadratic diver- 
gences and the finite contributions are scheme dependent 
and, in case of the if|-based calculational scheme, they 
are not gauge-invariant. Of course, this feature does not 
contradict the gauge-invariant renormaliz ability of the 
theory, but it makes the finite results scheme dependent. 
One has to remember that the effective action, in any 
particular dimension, is a sum of the series of the terms 
with flo: (i-i, 0-2, (23 etc. Therefore, we can conclude that 
the effective action is scheme dependent in any particular 
dimension. 

The scheme dependence which we have explored (mul- 
tiplicative anomaly), should be seen as a particular case 
of the nontrivial reaction of the effective action of quan- 
tum fields to the change of variables. The fact that such 
dependence can be relevant, has been known for a long 
time (see, e.g., [55]). It is interesting to see what is the 
consequence of this fact for the sum of the Schwinger- 
DeWitt series. We perform the corresponding calculation 
in the next section. 



IV. DERIVATION OF THE ONE-LOOP 
FORMFACTORS 

The effective action (21) can be expressed by means of 
the proper time integral, involving the heat-kernel K{s) 



sTr K{s). 



(50) 



of the operator O, 

fi=-i r ^ 

Let us derive this expression using the method developed 
previously in [56, 57] (see also [58] for the general review). 
According to Ref. [56] , we can perform calculations either 
by using Feynman diagrams or through the heat-kernel 
solution, which was obtained in [59, 60]. The result of the 
two methods are going to be the same, but the covariant 
heat-kernel solution is much more simple, so we will use 
this approach. 

The expression (50) can be expanded into the powers 

of the field strengths (curvatures), Rf^uap, S^u and P. 
Up to the second order in the curvatures the expansion 
has the form [59] 



Tr:K{s) = 



,4-2a 



I 



+ s''[R^^fi{-sD)R'''^ + Rf2{-sD)R 
+ Ph{-sn)R + Pf4-sn)P 
+ S^^M-sa)S''-^]}. (51) 



where the expressions for F and S^i, are defined in 
(34) and (35), u is the dimensional regularization 
parameter, fj, is an arbitrary renormalization parameter 
with the dimension of mass and the functions are 
given by 



Mr) 
/2(r) 

/3(t) 

h{r) 



fir) - 1 + r/6 

fir) ,^/(r)-l /(r)-l + T/6 



288 



+ 



24t 



/(T)^/(r)-l 



12 
l-/(r) 
2r 



2r 



8r2 



/4(r) = 



(52) 



where 



fir) 



-f 

Jo 



^^ga(l-a)r^ 



T = -SD . 



After some algebra, we arrive at the integral represen- 
tation for the effective action up to the second order in 
P, S^v and R^v, in dimensional regularization, 



10 



+ l^^ e'^RMk A^A'^ + IIR^, Mk R^" + lkRMkR], (53) 

where 

M - •^(*") M - -^(^"^ M - ■^(^^) M - ^ A/ - ^ 



and the nonzero coefficients for have the form 



k = - 



16 



1 

4 
1, 



'3 



/3 = - 

1, 



h = 



24 



7* 



1 

12 



/5 = 



7* 
'5 



2, T'^' 



-2. 



Indeed, the coefficients If^ and are the same as for the 
free fermion field [57], while the last set of coefficients is 
due to the interaction with A^ and do not have analogues 
in the free field case. For we have 





1 




1 


'4 




1 




4' 


h — ~ 


"2 ' 




2 ' 




]*DA 


1 


7* DA 

h — 




i*DA 




1, 




n — 


~2' 


-1 


1 I'll 






iDA 

h — 


-2, 


If^ = 2 












A2 = 


-1, 


A4 = l, 




A^ = -l, 




a: 


= 1, 


It"" = 


1 


jAR 

h — ~ 






1, 






~6' 


1, 


'4 — 








1 

6' 


iRA _ 1 




/I''' = -1 






(54) 



As far as the purely gravitational terms were calculated 
in [57], here we are mainly interested in the one-loop 

corrections to the term F^,^. Let us first perform all the 

calculations for the operator H^. Substitute (38) into 
(51), then the result into (50), and take into account 
only the terms with and f^. Using the notations [56], 



t = sM"^ , u = 



t 



□ 

'M2 



y = 1 - - In 
a 



/ 2 + a \ 



4n 



□ -4M2 

and the known results for the integrals [57], 

m2 \'^-2 r°° . , f{tu) 



(55) 



ri: 



dt e 



i-i)(^i)-i4.i 



/ M2 

a2 - 4 / I 



dt e 



3a2 
1 



18 



+ 0(2 - w) , 



+ 1 



4a2 



where we also denoted 
1 1 



+ 0(2-w). 



+ 2a) + 0{2 - w) , 



After some algebra, we arrive at the explicit expression 
for the one-loop correction to the classical F^^-term, 



d^x^F, 



- + kr{a) ^^^^(56) 



2(47r)2 

where the electromagnetic formfactor has the form [74] 



k 



r(a) = y(2 



3a' 



(57) 



A simple way to check our result is to compare it to 
the expressions for the Schwinger-DeWitt coefficients de- 
rived in the previous section. Let us remember that the 
calculation of the formfactor presented above has been 
performed in Ad. There is nothing that can be compared 
to ai , just because the quadratic divergences vanish when 
the dimensional regularization is used. The comparison 
to 02 is absolutely successful and complete, because here 
we meet two perfect correspondences. First, it is easy 
to see that the coefficient of the pole term of (56) is ex- 
actly the 02 from (45) at Ad. Second, it is easy to check 
that the UV limit a ^ of the formfactor (57) perfectly 
corresponds to the coefficient of the pole term of (56). 

Furthermore, we can partially compare our formfactor 
(57) to the expression for the 03, 1 in (48). Unfortunately, 
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the complete verification is not possible because certain 
part of (48) correspond to the next, third, order in the 
heat-kernel solution, which we do not use here. So, we 
can compare only the i^VVF-type terms. If we expand 
the corresponding part of fcf ^(a) into power series of the 
parameter a, in the IR limit, where 



a ~ 



□ 

M2 



use Maxwell equation 







(58) 



(59) 



and integrate by parts, we obtain the corresponding term 
in the effective action (48), 



which is actually well-known from [3]. 



(60) 



Consider the second calculational scheme, based on the 
operator H2, defined in Eq. (31). The calculations are 
a bit more involved, but are in general analogous to the 
previous case. Finally, we arrive at the following result 



^(1) 



2(4 



Y 



2 

+ 3 



3e 

8Y 



^ V3 3a2 



+ R 



Y 



V3a2 



R 
(61) 



1,FF 



(62) 



Once again, the UV test and the comparison with the 
03,2 are perfectly successful. However, the low-energy 
expression is not (60) anymore in this case, instead we 
have 



1607r2 



(63) 



As we will sec in Sect. 6, this divergence between (60) 
and (63) results also in the ambiguity in the decoupling 
theorem. 

It is possible to make an interesting verification of the 

results for the two formfactors fcf ^ (a) and fcf ^ (a) using 
the coefficients 0,3,1 and 0.3,2- If we expand (57) to the 
first order in operator □, we have 



k(^{a)=Y[2- 
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2 □ 



2 ^ 

3o2y 9~~15M2' 
that gives us a term in f which has the form 



f(i) 



15 (47r)2 



(64) 



(65) 



Using (59) in (65), we can see that 
2e2 



f(i) 



15 (47r)2 

^^^^^^y.^.V^V.F-V„F,«(66) 



and that is the same contribution given to f in M by 
the third term of Eq. (48), as it should be. To the form 
factor k2^{a), we find 



f(i) 



10 (4 



Now if we take the contribution of (49) to f in 4d the 
coefficient proportional to is the same as (67). 

It is important to note that the coefficients of the terms 
Ai^ Af, in Eqs. (66) and (67) are different in M. They 
are equal only in 6d as it was explained in Sect. 3. 



V. MASSLESS LIMIT AND TRACE ANOMALY 

In order to better understand the physical sense of 
the result (56) and also its relation with the MS-based 
anomaly-induced action (17), let us consider the high- 
energy limit. As was already explained in the previous 
sections, this can be done by taking the vanishing mass 
limit. On the other hand, this limit helps to establish the 
relation between the effective action (56) and the confor- 
mal anomaly (7). 

Consider the high-energy limit, when the mass of the 
quantum field (i.e. of electron) is negligible. Taking the 
limit o ^ 2 in the expression (56) with either one of 
the two available form factors, we arrive at the following 
leading-log behavior of the electromagnetic sector: 



(68) 



Similar asymptotic behavior takes place also in the gravi- 
tational sector of the theory. For instance, the Weyl term 
has similar form factor [56, 57], 



/3i C^^"/^ In (-^) 



(69) 



Let us note that the same asymptotic behavior can be 
recovered in the Minimal Subtraction - based scheme of 
renormalization [42, 61] (see also [62] for an alternative 
consideration), which is completely reliable in the mass- 
less case. 

The expressions (68) and (69) are sufficient to derive 
the corresponding parts of the conformal anomaly, even 
in case of a local conformal symmetry. For this end, let 
us apply the conformal parametrization of the metric (4) 
and the differential relation (2). Consider the case of (68) 
as an example, (69) is completely analogous. If we replace 
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the parametrization (4) into (68), the only place where 
the a field shows up is the □. This operator becomes 



□ = e" 



-2<T 



□ +0{d(j) 



(70) 



where we denote □ the d'Alembertian operator con- 
structed with the metric g^^, 



□ 



The explicit form of the terms 0{da) is in fact irrel- 
evant for us. When we apply (2), only the first term in 
the bracket (70) is important, because the other terms 
vanish after we set cr ^ 0. Of course, the logarithmic 
dependence makes 



In ^ = -2c7 + In — '- 

Finally, after applying (2) we arrive at 



(71) 



(72) 



in the electromagnetic sector. In a similar way one can 
obtain the wC^-term in the general formula for anomaly, 
Eq. (7). One has to note that the OR term and the 
Gauss-Bonnet term can be also derived from the local 
and nonlocal finite parts of effective action [59, 63]. At 
the same time, the OR part is a subject of an important 
ambiguity. The origin and mechanism for this ambiguity 
has been explained recently in [64] (see also [36] for the 
review). 



VI. RENORMALIZATION GROUP, 
LOW-ENERGY LIMIT AND DECOUPLING 

The two form factors (57), (62) contain all necessary 
information about the scale dependence of the coupling 
parameter e at the one-loop level, within the correspond- 
ing calculational schemes. In the present section we 
shall mainly restrict our attention to the flat space case 
and discuss this dependence in detail. Namely, our task 
here is to calculate the "physical" beta-functions in the 
momentum-subtraction renormalization scheme and look 
at their UV and IR limits. 

In the MS scheme the /^-function of the effective charge 
e is defined as 



de 4 e-^ 

/?e(MS) = lim u— = 

dfl 6 (47r)^ 



(73) 



The derivation of the /?-functions in the mass- 
dependent scheme has been described, e.g. in [65, 66]. 
Starting from the polarization operator, one has to sub- 
tract the counterterm at the momentum = , where 
M is the renormalization point. Then, the momentum- 
subtraction /3-function is defined as 



/3e 



hm M — 

n^4 dM 



(74) 



Mathematically, this is equivalent to taking the deriva- 
tive (we also write the same operation in terms of a) 



d , , ■2\C' d 
-ep— = e(4-a )- — 
dp ^ 'Ada 



(75) 



of the form factors in the polarization operator. If we 

apply this procedure to the form factor fcf^(a) of the 
f^j^-term, the expression for the /3-function in a mass- 
dependent scheme is 



47r)2 I 



6a3 (47r) 



20a^ - 48a 
2- 



(76) 



that is the general result for the one-loop /3-function valid 
at any scale [75]. 

As the special cases we meet the UV limit p^ » m?, 
or a ^ 2, 



1 UV 



4p/ 



3 (47r)2 



O 



(77) 



that is nothing else but the MS scheme result (73) plus 
small correction. In the IR regime, however, when p^ <C 
m^, the result is quite different, and moreover depends 
on the calculational scheme. For the first case H^, we 
have 



1 IR 



(47r)2 



4M2 
15 m2 



+ o(^). 



(78) 



This is exactly the standard form of the decoupling the- 
orem [32]. 

Similar calculations starting from the form factor 
k^^ia) give 



= 12^ + 
-3(a^-16)ln(|^)}. 



(79) 



In the UV limit p^ 3> m? , the above /3-function is in 
agreement with the standard result (77), while in the IR 
limit p^ <^ we obtain 



0, 



2 IR 



5(47r) 



(80) 



As we can see from the expressions (78) and (80), there 
is a slight difference in how /3-functions go to zero in the 
IR limit. To discuss the physical sense of this fact, let 
us take the difference between the two form factors (57) 
and (62) 

AF^^ = fcr(a) - kria) = a[i - ±) - ^ (81) 

and expand AF^^ in power series of a. In the IR limit 
we have ^ p^ / M"^ = - D/M^, so 



AF^^ = 



1 



□ 



30 M2 



0{p^/M^) . 



(82) 
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This is exactly the difference in the form factors which 
caused the ambiguity (calculational scheme dependence) 
in the decouphng theorem. In order to evaluate the 
source of this difference in the effective action, we should 
consider what would be the new terms in the equation of 
motion, generated by the term 

As this term is proportional to the operator □, if we are 
working in flat spaces, we use Eq. (59) to obtain a term 
proportional to (V^F^")^. This term will not influence 
the equations of motion in flat space in the O(e^) ap- 
proximation [3]. However, as we will discuss later on, the 
situation can be different in curved space. 

VII. RUNNING CHARGE 




0.0 05 1.0 1.5 , 2.0 2.5 



FIG. 1: The beta functions corresponding to form factors fci 
(solid line) and k2 (dashed line). The vertical axis is drawn 
in units of e^/(47r)^. 



Let us now consider the running of the electromag- 
netic charge. The difference in the expressions for the 
form factors k[^{a) and /c|"^(a) and eventually in the 
expressions for the /3- functions (76) and (79), means that 
the effective charge (i.e., the running charge) may have 
different behavior for the different calculational schemes 
(that means, for the ones based on and H2 operators). 
Also, we expect that the running within the Momentum 
Subtraction scheme will be distinct from the one within 
the Minimal Subtraction scheme, especially in the low- 
energy region. Let us check out what the real situation 
is. 

To investigate the renormalization of the correspond- 
ing quantities, in the physical scheme, let us apply the op- 
erator ~epd/dp to fcf^(a) with □ traded for — p^. In do- 
ing so, we find the expressions for the /3-functions, which 
can be conveniently presented as 

Here t is a dimensionless parameter defined hy t = 
In {p/m). 

The UV limit is achieved for p 3> m, or equivalently 
t S> 1, while in the IR limit we meet inverse relations 
p <^ m and t <^ 1. With the help of the Mathemat- 
ica computer software, one can calculate explicitly the 
above beta functions and integrate them. For a clear il- 
lustration, we plot the beta functions for both cases as 
fimctions of the parameter a (see Fig. 1). 

The integration of the renormalization equation cor- 
responding to /3i can be performed by using the Math- 
ematica software such that the integral curve describes 
the running coupling constant in the physical scheme. 
For comparison, we plot in Fig. 2 this curve together 
with the curve for the running parameter for the Minimal 




FIG. 2: The curves for the running parameter e(f)^ seem to 
coincide for both MS and physical scheme (in the last one we 
use the prescription corresponding to the form factor fci). Ac- 
tually, both curves are slightly different as well as the Landau 
pole (see Fig. 3). We have used e(0) — 0.1. 



Subtraction scheme, for large value of t, where a Landau 
pole shows up. It is not easy to visualize a difference 
between both cases, which look almost identical. Actu- 
ally, the curve for the running parameter in the physical 
scheme is shifted a little bit to the right. This difference 
can be made clear if we increase the plot scale in the re- 
gion around t ^ 5900, as illustrated in Fig. 3. Here, one 
can see actually two vertical lines (the dashed one for 
the MS scheme), indicating two different Landau poles. 
They are, however, very close, since the corresponding 
values for t differ by less than 0.02%. The situation is 
very similar to the one described earlier for the scalar 
field theory [68]. 



The running of the effective electromagnetic charge is 
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FIG. 3; The same plot as in Fig. 2, viewed with more detail 
around t ~ 5900. Here the two vertical lines indicate different 
Landau poles. The dashed line corresponds to the MS scheme. 
The value for t corresponding to the poles differs from each 
other by less than 0.02%. 



FIG. 4: Curves for e{t)~ . The dashed line corresponds to 
the MS scheme and the continuous line for the momentum 
subtraction scheme. For higher t, both curves are parallel 
straight lines. Substantial difference between the two plots 
takes place for t < 1. 



shown in Fig. 4, where we plot versus t = In {p/m). 
The plot for the MS case is a straight line, as it has to 
be. The plot for the momentum subtraction scheme rep- 
resents two straight lines in the asymptotic (IR and UV) 
limits with the smooth transition between them in the 
intermediate region. It is easy to note that the curves 
are very similar to each other in the UV region, that is 
for t ^ 1 (in fact this feature holds already for t > 2). 
At the same time, even in the UV the two plots do not 
coincide and are represented by two parallel straight lines 
with slightly different initial points (at t = 0, that means 
p = ^0 and = HQ. The effect of the effective UV shift 
of initial point will take place not only in QED, but also 
in QCD and EW sectors of the SM and beyond. This 
effect, despite being very small, may have interesting ap- 
plications. For example, one can take into account the 
effective shift of the initial values of the couplings when 
calculating the UV effects in gauge theories via the renor- 
malization group (see, e.g., [69]). The calculations of this 
sort are quite relevant, including for the perspective LHC 
physics. Taking the effect of the masses of the quantum 
fields into account may, in principle, improve the preci- 
sion of the results. As a particular case, the effect of the 
initial data shift can be also observed in the supersym- 
metric models such as MSSM and its extensions. As a 
result, there may be a slight violation of the exact con- 
vergence of the running coupling constants 5, g' and gs 
which will likely form a small triangle rather than meet 
in a single point. This feature of the massive theories has 
been already discussed in [70] , but the effect can be done 
more explicitly when using analytic expressions for the 
/3-functions. 



VIII. ON THE CONFORMAL VIOLATION IN 
MASSIVE CASE 

In the previous sections we have considered the exist- 
ing interface between the quantum contributions coming 
from the massive and massless field loops. In particular, 
we have seen that, in the high energy limit, the form fac- 
tors and the /3-functions tend to the ones in the massless 
case, that means they become close to the ones in the 
minimal subtraction scheme of renormalization (MS, in 
our case). As we already learned in section 2, in this 
situation the violation of local conformal symmetry oc- 
curs due to anomaly and the corresponding EA can be 
presented in a closed form (17). 

On the contrary, the low-energy quantum effects are 
characterized by the phenomenon of decoupling [32] . As 
a result, in the IR limit the quantum effects are quadrati- 
cally suppressed and one can rely on the classical Maxwell 
equations. An interesting question is what is the differ- 
ence between the violation of local conformal symmetry 
in the "far UV" limit and in the "far IR" limit, where we 
have a remnant quantum effects. Let us emphasize that 
there is no contradiction in having two different forms 
for the violation of conformal symmetry, because they 
correspond to the two different physical situations. 

Let us first compare the dependences on the confor- 
mal factor in the two symmetry-violating terms. For the 
sake of simplicity we can take a — const, that means 
the global conformal symmetry. The situation for the 
local case will not be very different. In the UV limit, 
exactly as in the purely massless case we have, according 
to the Eq. (71), the linear dependence on a. Of course, 
the same result follows from the complete expression (17) 
derived for the precisely-massless case. On the contrary, 
in the situation of IR decoupling, we can take the low- 
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est order terms directly from the as-coefficient. In this 
case one meets the FDi^-type expression (83) and, also, 
RFF-typc terms. In both cases the cr-dependence is ex- 
ponential, that means the symmetry is violated by the 
terms which have a scaling law e~^'^. So, it looks like the 
scaling of the symmetry-violating terms is even stronger 
for the low-energy sector. However, this is nothing but a 
wrong impression. 

One can view the situation from a different position, 
if analyzing this question using the calculated expres- 
sion (56). Then the physical sense of decoupling in the 
low-energy limit becomes much more explicit. If we com- 
pare the FDi^-type expression (83) and the classical term 
F^j^, it is clear that the former has an extra factor of 

jm^, where is the square of the momentum of the 
photon and the rrfi is the square of the mass of the elec- 
tron. In the low-energy limit <C m? , hence we meet 
the simplest (and very clear) form of the Appelquist and 
Carazzone [32] quadratic decoupling law. In the case of 
RFF-type terms, for most of the physical situations, the 
decoupling is even much stronger. The reason is that 
the scalar curvature is proportional to the square of the 
typical energy of the gravitons and this energy is much 
smaller than the one of the photon. For instance, in the 
cosmological setting we have R oc H'^, where H is the 
Hubble parameter. In the present-day Universe the cor- 
responding values are Hq « 10"'*^ GeV for the Hubble 
parameter and e = -\/^ ~ 10~^^ GeV for the energy of 
the CMB photon. Similar relation between the two quan- 
tities holds during most of the evolution of the Universe. 
In this situation, the most important low-energy con- 
tribution to the violation of conformal symmetry comes 
from the (83) term, which has a nontrivial flat limit. 

IX. DISCUSSIONS AND SUMMARY 

The two distinct approaches to the derivation of one- 
loop quantum corrections to the photon sector of the 
curved-space QED have been explored. First, we derived 
the anomaly-induced action, coming from the integration 
of conformal anomaly. The effective action which we gain 
in this way has an enormous advantage of being exact, 
in the sense it is not related to some particular form of 
series expansion, except the one into the loops (taken in 
the first order). The representations obtained here can 
serve, in principle, as a consistent background for the in- 
vestigation of quantum processes in the early universe, 
when the masses of quantum matter (fermionic) fields 
are irrelevant. The anomaly- induced action has yet an- 
other advantage of being scheme independent, because 
it is based on the Minimal Subtraction renormalization 
scheme and is, after all, controlled by the one-loop diver- 
gences. 

On the other hand, in the massive case there is the 
violation of conformal symmetry coming from the form 

factor in the electromagnetic sector of QED in curved 
space-time. This calculation is based on the physical 



renormalization scheme and hence it is supposed to be 
more adequate in the later universe, when the masses of 
the fields play an important role. The price one has to 
pay for a more consistent physical approach is related to 
the restricted power of the available calculational meth- 
ods, which are equivalent to the use of common Feynman 
diagrams for the linearized metric perturbations on the 
flat space-time background. The corresponding result 
which we obtain here includes terms which are quadratic 
in F^j^ and may also depend on the curvature tensor. It is 
complete in the case of a flat space-time background, but 
it is not supposed to be a complete one for the curved 
space-time, where we can expect many higher order in 
curvature corrections which can not be calculated ex- 
actly. An interesting point is that we have found that 
the quantum correction depends on the choice of the cal- 
culational scheme [33]. Thus we have proven the exis- 
tence of the nonlocal and renormalization independent 
MA in quantum field theory. One of the consequences 
of this anomaly is the ambiguity in the prediction of the 
Appelquist and Carazzone theorem [32], which provides 
two different coefficients of the quadratic decoupling law 
at low energies. 

Since the quantum corrections in the electromagnetic 
sector include some ambiguity in the IR region, one 
should ask which one of the two schemes gives a cor- 
rect result. In our opinion the advantage should be given 
to the one derived through the operator H*, because it 
is more natural and preserves gauge invariance. How- 
ever, it is worthwhile to be aware of the ambiguity which 
is a manifestation of a typical property for the off-shell 
effective action in quantum field theory. In the present 
case this ambiguity becomes essential due to the pres- 
ence of an external gravitational field. One can note that 
the simpler form of quantum correction (17), derived via 
conformal anomaly, is also ambiguous due to the pres- 
ence of an arbitrary functional Sc- Of course, the two 
kind of ambiguities are unrelated, but they can be seen 
as manifestations of a general feature of effective action. 
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Note Added. 

After this paper was resubmitted, we learned about the 
recent paper [71] (see also [72]) where the same problem 

was considered using technically different approaches. 
We would like to cite this paper, especially because the 
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parts of conformal anomaly-induced action arc different 
only due to the choice of paramctrization for the auxil- 
iary fields. Also we would like to point out that, despite 
of the mentioned differences, the results are qualitatively 



similar. In particular, in both approaches the conformal 
anomaly can be partially restored from the calculations 
of massive loop in the massless limit. 
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